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Abstract 

An original boundary integral formulation is proposed for the problem of a semi-infinite 
crack at the interface between two dissimilar elastic materials in the presence of heat flows 
and mass diffusion. Symmetric and skew-symmetric weight function matrices are used to¬ 
gether with a generalized Betti’s reciprocity theorem in order to derive a system of integral 
equations that relate the applied loading, the temperature and mass concentration fields, the 
heat and mass fluxes on the fracture surfaces and the resulting crack opening. The obtained 
integral identities can have many relevant applications, such as for the modelling of crack 
and damage processes at the interface between different components in electrochemical en¬ 
ergy devices characterized by multi-layered structures (solid oxide fuel cells and lithium ions 
batteries). 
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1 Introduction 


Modelling of interfacial crack problems in elastic bimaterials in presence of thermodiffusion rep¬ 
resents an important issue for many engineering applications. In particular, it is crucial for 
studying fracture initiation and propagation at the interface between different components of 
electrochemical energy devices which are characterized by multi-layered structures, such as solid 
oxide fuel cells and lithium ions batteries. Indeed, due to the high operational temperature that 
can be reached and to the intense particle fluxes that are required for maintaining the electrical 
current, the components of these devices ar e subject to severe thermomechanical stresses as well 


as stresses induced by the particle diffusion ( Anandakumar et ah . l2O10l : [^ang et ai] . l2O10l i. This 


can cause damage and crack formation co mpromising their performan c es in terms of pow e r gen¬ 


eration and energy conversion efficiency ( Lowrie and Rawlings . 2000; Malzbender et ah . 20031 : 


Goutianos et al 


201fll : Zhao et ah . 2010; Pharr et ah . 2013l l. For this reason, the modelling o:’ 


fracture processes at the interface between the components of such battery devices is fundamen¬ 
tal to the prediction of these phenomena, and subsequently enables successful manufacture and 
reliable performances of the systems. 

Due to their distinctive feature of reducing by one the dimension of the considered prob¬ 
lems, so that only the boundary or surface of the domain needs to be modeled, boundary 
integral formulations are particularly suitable for studying multiphysics phenomena such as 
dynamic and static fracture processes in thermodiffusive elastic materials. Several bound¬ 
ary integral approaches have been proposed for solving crack problems in linear elastic, ther¬ 
moelastic and thermodiffusive materia l s, suc h as analytical techn i ques b ased on the metho d 


of singular integral equations (jWeaveii . 119771 : iBudianskv and Ricel . Il979l: iLinkov et al.l. Il997li 


and numerica l techniques based on the boundary e lement method (jRizzo and ShipovL [1977 


Brebbia et ah. 

00 

Sladek and Sladek. 

1983, 

1984a 


1984a 3)- In ah these approaches, the displace¬ 


ments and stress fields are defined by integral rel ations involving the Green ’s functions, which 
need to b e derived analytically in explicit form (|Bigoni and Capuanil . l2002l l or computed nu¬ 
merically ( Ang and Telles . 20041 1. Altough Green’s functions have beei i derived for several crack 
problems in linear thermoelastic and thermo diffusive elastic materials ( Sturla and Barber . 19881 : 
Hou et al. , 20081 : Kumar and Chawla , 2012a 3) , their utilization for calculating physical displace¬ 


ments and stress fields on the crack faces requires challenging numerical estimation of integrals 
whose convergence should be assessed carefully. Moreover, the approach based on the Green’s 
function method works when the tractions and the thermal and diffusive stresses acting on 
the discontinuity surface are symmetric, but not in the case where asymmetric mechanical and 
thermodiffusive loading distributions are applied on the crack faces. 

In this paper, the problem of a semi-infinite quasi-static crack at the interface between two 
dissimilar thermodiffusive elastic materials is addressed by means of an origianl boundary inte¬ 
gral formulation which avoids the use of the Green’s functions and the challenging computations 
connected, without any assumptions regarding the symmetry of the loading and of the tempera- 
ture and mass concentra t ion p r ofile at the interface. The general approach r ecen tly proposed in 
Piccolroaz and Mishurid ( 2013l k Morini et al. ( 2013al k Vellender et al. ( 2013l l and Mishuris et al 


( 20141 ') for inter facial crack problems in isotropic and anisotropic elastic bimaterials, based on 
Betti’s reciprocal theorem and weight functions theory, is extended in order to study fracture 
processes in presence of thermodiffusion. The volume integral terms present in the reciprocity 
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identity, associated with the temperature and mass concentration effects ( Nowacki . 1974al i. are 
converted into su rface integrals th rough an exact transformation based on the notion of Lame 
elastic potentials I Slaughter 120011 1 while assuming that the temperature and mass concentration 


are harmonic in the domain. The derived original form of Betti’s ident hy is used together with 
symmetric and skew-symmetric weight function matrices derived by Piccolroaz and MishurisI 
( 20 Til l for formulating the considered crack problem in thermodiffusive bimaterials in terms of 
singular integral equations. 

The article is organized as follows: in Section [2l the static governing equations for a lin¬ 
ear elastic thermodiffusive media are formulated. Reciprocity indetities are introduced and the 
volume integral terms associated with the temperature and mass concentration fields are trans¬ 
formed into boundary surface integrals by introducing Lame elastic potentials and applying 
second Green’s theorem. In Section [3l the problem of a quasi-static crack at the interface bew- 
teen two dissimilar elastic thermodiffusive media is introduced. The weight functions, defined 
as a special singular solution of the homogeneous traction-free problem are used together with 
the obtained Betti’s identity for formulating the problem in terms of boundary integral equa¬ 
tions. The case of a plane strain crack is analysed in Section SI Lame potentials in the weight 
functions space are derived in closed form, and explicit integral identities relating the applied 
mechanical loading, the profiles of temperature, mass concentration, heat and mass fluxes on 
the interface and the resulting crack opening are obtained. Finally, in Section [5l the integral 
identities are used to study some illustrative examples of plane crack problems in the presence 
of thermodiffusion. Exact expressions for the crack opening and tractions ahead of the crack 
tip associated with the introduced temperature and heat flux distributions are derived, and the 
corresponding stress intensity factors are calculated in closed form. 


2 Preliminary results: governing equations and reciprocity the¬ 
orem 

In this Section, constitutive relations and static balance equation for linear thermodiffusive 
elastic media in infinitesimal deformations are introduced. Betti’s integral identities are derived 
by means of the theorem of the reciprocity of work, and volume integral terms associated with 
heat flows and mass diffusion are converted into surface integrals using an exact transformation 
based on the introduction of Lame elastic potentials. 


2.1 Governing equations 

For a linear isotropic elastic body where temperature changes and mass diffusi on are considered , 
the constitutive relationship between the stress <t and the strain e is given by ( Nowacki . 1974al l: 


<T = 2/j.s + (Mrs - jtO - 7cX) I, 


( 1 ) 


where A and /r are Lame’s constants, 9 = T — Tq is the temperature of the medium with respect 
to the reference state (e = 0, T = Tq C = Cq), x = C — Cq \s the concentration of diffusing 
particles with respect to the natural state, / is the identity matrix, xt = (3A -|- 2fj.)at and 
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7c = (3A + 2fi)ac with at and ac coefficients of thermal and diffusive expansions, respectively. 
Then in the natural state of the system we have 


e = 0, 0 = 0, x = 0, (2) 

and in order to desc ribe the material using linear thermoelastic diffusion theory ( Nowackl . 1974al : 


Sherief et al.l . 12004 ') the values of 0 and x are assumed such that \0/Tq\ <C 1 and Ix/CqI <C 1. 


In the static case, the equations of equilibrium are given by 


V • (T + 6 = 0, (3) 

where b represents the body forces. Substituting the constitutive relation ([I|) in ([3]), the Navier’s 
equations in terms of displacements u are obtained 

+ (A +/r)V(V • ri) + & = ytVS + 7cVx. (4) 

The conservation of energy and mass in steady-state conditions leads to 

V ■ q = W, V ■ j = 0, (5) 

where W is the heat generated for unit of time and volume by internal sources. For the considered 
linear isotropic thermodiffusive media, the heat flux q and particles current j can be expressed 
as 

q = -hVO, j = -DcVx, (6) 

where kt is the thermal conductivity coefficient and Dc is the mass diffusivity of the material. 
Using relations ([6]) in ([5]), the following Poisson’s and Laplace’s equations for the steady-state 
temperature and mass concentration fields are derived: 

ktM + W = 0, Ax = 0. (7) 

2.2 Reciprocity theorem: transformation of the volnme integrals 

Let us consider a thermodiffusive elastic body of volume V subject to the action of body forces 
tractions heat sources surface heating to the temperature 0^^^ and mass concen¬ 

tration on the boundary x[^^ • These causes are written symbolically in the compact form 

xa) = |&(i),t(i),iTa),0W,;^w}, (8) 

and produce in the body the state characterized by displacements temperature 0^^\x) 

and mass concentration 

pT) = 0^^\x), x*'^H®)| j x ^ V. (9) 

The stresses and the strains are connected by the relationship ([1]), where the dilatation 
is defined as tre = V • ri, and they are assumed to be continuous together with their first 
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derivatives. The displacements, temperature and mass concentration are also continuous and 
moreover have continuous derivatives up to the second order for x + dV. These functions 
satisfy the fields equations (l3|), ©(i) and (ED (2) with boundary conditions: 

t^^\x) = a^^\x)n{x), X £ dV, (10) 


and 


e^^^x) = el^\x), x£dv, 


( 1 )/ 


( 11 ) 


where n denotes the unit outward normal to the surface dV. 
Introducing another set of causes and effects 


IP) = .xf} , CPI = x'=>(x)} , (12) 


and applying the procedure illustrated in Nowacki ( 1974al Ol for elastic thermodiffusive media, 
the following reciprocity integral relations bewteen the two systems of causes and results are 
derived; 


/ 

Jav 


cr(^^n • ■ u 


+ j -^(2)-6(2) .1^(1)^ 


+ 


It J -0(2)v.uW)dl/ + 7,y -;^(2)v.iu(i))dl/ = 0, (13) 


lav 


(Vx(i)y(2) _ vx(2)^(i) 


n dS = 0. 


(15) 


Expressions ()13p and (1141) have been extensively used in the literature in order to develop 
numerical buondary elements methods for s tudying both static and dynamic c rack problems in 
elast ic materials subject to thermal stresses ( Sladek and Sladek . 1983I . 1984a b; Dell’Erba et al. 


1998l l. In most of these approaches, starting from equations (fTH|) and (fTTl) . integral expressions for 
the physical displacements, stresses and temperature and are derived as functions 

of the test quantities <t(^) and which in general are assumed to be fundamental solution s 


of the evolution equations or weight functions ( Brebbia et ah . 1984 : Rizzo and Shipov . 1977l l. 


Several analytical and numerical transformations have been proposed for reducing the volume 
integrals associated with body forces and coupli ng between mechanical strains and temperature 
to surface integrals ( Shiah and Tan . 19991 . 2012 1. 

Here, we assume zero body forces and heat sources acting on the system, and we introduce 
an exact procedure for transforming coupling volume terms involving temperature and mass 
concentration in equation (I13|) . Erom Helmholtz’s decomposition theorem, since both the dis¬ 
placements are solutions of the equilibrium equation ([3|), they can be repr esented by 


scala r displacement potentials and vector displacement potentials (Slaughte 

l 2 nnil l: 

-g V X u(2) = Vg)(2) -g V X 0(2), (16) 
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substituting expressions (fTHjl into CSD, the following equation is obtained: 


lav 


) dS+ 


'yt J dF + 7c y dF = 0. (17) 

The Green’s second identity states that two arbitrary scalar functions, (p and id, must satisfy 

[ {(pM - ^A(P) dV = [ {(pV^ - i9V(p) ■ n dS. (18) 

Jv Jav 

Substituting respectively to and to (p, using equations dZD and 

remembering that zero heating sources are assumed, we get 


iv 


iv 


lav 


lav 


,(2)_^(2)v0(l)y 

■ n dS, 

(19) 

,(i) _ (^(i)v6»(2)^ 

■ n dS, 

(20) 

>< 

> 

1 

• n dS, 

(21) 

,(1) _ 

■ n dS, 

(22) 


Jv Jav ^ 

substituting expressions (fT^ - ([2^ into (HD, we finally derive 

f (cr^^^n ■ — cr^‘^^n ■ dS+ 

+7t j ■ n dS+ 

+ 7c / _ (^(2)vx(i) - x(2)V(^(i) + .ndS = 0 . (23) 

Jav ^ ^ 

By means of the proposed general procedure, the volume integral terms associated with ther¬ 
mal and diffusive stresses in the reciprocity identity (fT7t) has been reduced to surface integrals. 
Remembering the fluxes definitions ([6]), expression (12311 can be written as follows 

f (cr^^^n-u^‘^'l—a^‘^^n-u^^'^'\dS+ 

Jav ^ 2 

-|-7t J ■ n dS + Pt J • n d^-l- 

+ 7c / ■ n dS + (3c [ ■ n dS = 0, (24) 

Jav ^ ^ Jav ^ 2 
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where j3t = 'yt/kt and /3c = JtlDc- Assuming W = 0, reciprocity identities (fH)l and (fT^ can 
also be expressed in terms of fluxes: 


/dv 


^(2)0(1) _q(i)0(i)^ ■ndS = 0, 


lav 






n dS = 0. 


(25) 


In the cases where it is possible to determine explicitly the elastic potentials (USD , us¬ 
ing expression (1241) together with (I25D boundary integral formulation of static crack problems 
i n thermodiffusive solids can be obtained avoiding numerical estimation of volume integrals 
( Cheng et alJ . 2001 b In the next Sections, the obtained relation (Oip will be extensively applied 
in order to derive integral identities for the modelling of fracture phenomena at the interface 
between dissimilar elastic materials in presence of thermodiffusion. 



Figure 1: Geometry of the model. 


3 Interfacial cracks in thermodiffusive media 


A static semi-infinite crack between two dissimilar elastic materials in presence of thermodif¬ 
fusion is condidered, the geometry of the system is shown in Fig(T] No body forces, heat and mass 
sources are assumed. The crack is situated in the half-plane = {jc = (xi, X 2 , X 3 ) G : xi < 0, X 2 = O}, 
and general non-symmetric loading applied to the crack faces is assumed. Further in the text, 
we will use the superscripts and “ to denote the quantities related to the upper and the lower 
elastic thermodiffusive half-spaces, respectively. The applied loading c an be decomposed in th e 


symmetrical and skew-symmetrical parts, which are defined as follows ( Piccolroaz et ah . 2009l l 


(P) = 2 + P ) ’ [[P]] =P^ -P > 


(26) 


where we used standard notations to denote the average, (/), and the jump, [[/]], of a function 
/ across the plane containing the crack, X 2 = 0, 

if) {xi.X'i) = ^[/(xi, 0 +,X 3 )-h/(xi, 0 ",X 3 )], [[fl{xi,X 3 ) = /(xi,0+,X 3 )-/(xi,0“,X 3 ). (27) 
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In the proposed model, thermal conduction and diffusion are both assumed to be isotropic, 
and then only normal stresses are induced by temperature changes and mass diffusion (see 
constitutive relation ([T])). 

Further in the paper, using the Betti identities (l2^ and ([25]), the considered crack prob¬ 
lem will be formulated in terms of boundary integral equations. In order to correctly apply 
expressions (I24p and ()25p . a preliminary discussion regarding the behaviour of the temperature 
and mass concentration functions in the upper and lower elastic thermodiffusive half-spaces is 
reported in next Section 13.11 and boundary conditions for these fields and the fluxes dH) are 
introduced. 


3.1 Temperature and mass concentration in thermodiffusive half-spaces 

Since zero heat and mass sources are assumed, the temperature 0 and mass concentration fields 
are harmonic functions determined by the solution of a Laplace’s equation assuming the form 
of d?]) ( 2 ) • Considering the upper elastic thermodiffusive half-space these equations become 

A0+ = 0, Ax+ = 0, X2> 0. (28) 


Note that the conditions of harmonicity of the temperature and concentration fields (1281) must 
be satisfied in order to transform the volume integrals involved in the reciprocity identity dni) 
in surface integrals following the procedure explained in Section 12.21 

The values of 6 and x for X 2 = O"'', corresponding to the plane containing both the bounded 
interface and the crack, are given by: 

0+(x2 = 0+) = 0(xi,O+,X3), x'^(a;2 = 0+) = x(a:i,0+,X3). (29) 


Replacing 0+ with 6~ and x^ with x~ a couple of Laplace’s equations identical to the ([28]) and 
defined for X 2 < 0 are obtained for the lower thermodiffusive half-space, and then for X 2 = 0 “ 
boundary conditions analgous to ([ 2 ^ can be written. 

The average and the jump of the temperature and concentration are defined on the whole 
plane X 2 = 0 : 


(0) (xi,X3) = i[ 0 (xi,O+,X 3 ) -F 0(a;i,O ,^3)], 

— 00 < Xl, 3:3 < -boo. 


{x) ( 3 ^ 1 , 3 : 3 ) = ^[x(3:i,0+,X3) -Fx(3;i,0“,X3)], 

— 00 < Xl, X3 < -boo. 

(30) 

[[0]](3:i,a:3) = 0 (xi,O+,X 3 ) - 0 ( 3 : 1 ,0“, 3:3), 

— 00 < Xl, X3 < -boo. 


[[x]](3:i,X3) = x(3:i, 0+,X3) - x( 3 :i,0",3:3), 

—00 < Xi,X3 < -boo. 

(31) 


Further in the text, we will refer to (13QD and dni) as to imperfect thermodiffusive interface 
conditions. The values of the components of the fluxes normal to the interface plane, q 2 = Q ■ ^2 
and j 2 = J • 62 , for X 2 = O'*' are 


q+{x2 = 0 +) = - kt 


de+ 


dx2 


dx' 


= 92(3^1, 0 +,X 3 ), j+(x 2 = 0+) = - I?;. 

0:2=0+ 


0:2=0+ 


= J 2 (a;i, 0 +,X 3 ), 

(32) 






and similarly to the cases of temperature and mass concentration, expressions analogous to the 
(j 32 p can be defined for the normal fluxes in the limit X2 = 0 “. 

The average and the jump of of q2 and j2 are also defined on the whole plane X2 = 0 ; 

(^2) {xi,x^) = ^[92(3:1,0+,X3) + g2(a;i,0“,a;3)], -00 < xi,x^ < +00, 

O2) (a;i,a;3) = ^b2(a;i,0+,X3) + j2(a;i,0“,X3)], -00 < xi,X3 <+00, ( 33 ) 

[[g2]](a;i,a;3) = g'2(a;i, 0 +,a;3) - ^2(3:1, 0 “, 3:3), -00 < a;i,X3 < +00, 

[[j 2 ]]( 3 :i,a: 3 ) = J 2 ( 3 :i, 0 +,X 3 ) - j 2 (a;i, 0 “,a: 3 ), -00 < xi, 3:3 <+00. ( 34 ) 

Both in the upper and lower half spaces, the temperatures and mass concentration fields are 
assumed to vanish at the infinity: 

6 ~^ = = 0 ) for xi^X'i ±00, 0 ^ = = 0 ,for X2 —)• +00, 


6 =x =0, for xi, X 3 —>■ ± 00 , 9 =x = 0,for X 2 —>■ — 00 . (35) 

Since the heat and mass conduction problems (1281) are stationary, the flux functions <72 
and j 2 must satisfy the self-balance conditions on the boundary. For the considered upper 
thermodiffusive half-plane this means that: 


-1-00 


r |■+Qa 

92(3:1, 0 +,X 3 )dxidx 3 = 0, / / 92(2:1,0+,X3)dxidx3 = 0, 

J J —00 


(36) 


where the definitions (I32p for the values of the fluxes on the interface plane x = 0 have been 
utilized. Note that the self-balance conditions for the fluxes on X 2 = 0“ are given by relations 
similar to the (|36p . Remembering the definitions (I33p and (I34p . integral balance conditions 
analogues to the (i36P can be derived for the average and jump of the fluxes across the plane 
X2 = 0. 


3.2 Boundary integral equations 

Considering the geometry of the model shown in Fig. [T| the Betti identities (I24p and ()25p are 
applied to a semi-spherical domain of radius r in the upper and in the lower half spaces and 
In the limit r —)• 00 , assuming that both the displacement fields and decay suitably 
fast at infinity and remembering conditions (1351) for the temperature and mass concentration 
fields, the reciprocity relation ()24p in the upper half-space becomes: 


'(0:2=0+) 


4^^(xi, 0+,X3) • n(^)(xi,0+,X3) - <Tf^(xi,0+,X3) • u’'^\xi,0'^,X3) 


+7r 


'(X2=0+) 


0 X 2 0 X 2 


+Pt 


I(X2=0+) 


(7?(^)(xi, 0 +, X3)4^^(xi, 0 +, X3) - (xi, 0 +, X3)4^^(xi, 0 +, X3) 


dxidx3+ 

dxidx3+ 
dxidx3+ 
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+7^ 


+ /^c 


'(X2=0+) 


X 




dxidx'i+ 


'(X2=0+) 


(xi, 0+, (xi, 0+, X3) - (xi, 0+, X3)jf ^ (xi, 0+, X3) 


dxidxs = 0 , 


(37) 

where u = [ui, n 2 , ^ 3 ]^ and <72 denotes the traction vector acting on the plane X 2 = 0 : <72 = < 762 - 
In the same limit r —)• 00 , expressions ()25p assume the following form: 


'(X2=0+) 


6'(^)(xi,0+,X3)4^^(a;i,0+,X3) - 6»(^)(a;i,0+,X3)g^^^(xi,0+,X3) dxidxs = 0, (38) 


[ X^^\xi,0+,X3)j[^\xi,0+,X3) - x^'^\xi,0+,X3)j[^\xi,0+,X3) dxidx 3 = 0. (39) 

J {X2=0+) ^ 

Considering the boundary a; = 0“ instead of a; = 0+, three similar expressions are de¬ 
rived for the lower half-space. In the cases where elastic potentials can be derived explicitly, 
integral equations (f37|) . (l38l) and (l3^ can be used together with their analogues in the lower 
half-plane in order to obtain expressions for the physical fields associated with the crack prob¬ 
lem (72^\ 0*^^) and These physical quantities are evaluated by means of a 

set of predetermined auxiliary functions a^ 2 \ Q^\ satisfying the balance 

equations. 


3.3 Betti formula and weight functions 

Following the approach proposed in sever al boundary element formulations of thermoelasticity 
( Chener et ah . 2001 : Shiah and Tan . 2012l b for the auxiliary solutions system we assume 0(2) = 
< 72 ^^ = 0 and = 0. The displacement field ri(2) is represented by the non-trivial 

singular s olution of t he homogen eous traction-free crack problem, commonly known as the weight 
function ( Bueckner . 1985, 1989l b and defined as follows 


u 


( 2 ) 


(X1,X2,X3) = RU{-Xi,X2, -X3), 


(40) 


where R is the rotation matrix: 


R = 


-10 0 
0 1 0 
0 0-1 


The transformation (1401) corresponds to a change of coordinates consisting in a rotation of 
an angle vr around the X 2 —axis, and it is connected to the fact that the weight function U is 
defined in a different domain with respect t o physical displacement, wh ere the crack is placed 
along the positive semi-plane xi = 0, X 2 > 0 (|Willis and Movchanl . Il995l i. The stress tension U 
associated with the displacements U is introduced: 


( 7 ( 2 )(xi,X 2 ,X 3 ) = RI:(-Xi,X2, -X 3 )R. 


(41) 
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Remembering the general definition (I16h and taking into account the transformation (I40p . the 
following elastic potentials are introduced in the space of the weight functions: 

ip^‘^\xi,X2,Xz) = ^{-Xi,X2, -X3), 't/^‘^^^(xi,X2,X3) = R^(-Xi,X2, -X3). (42) 

Substituting these expressions into the Betti formula (1571) . and replacing U{xi,X 2 ,X 3 ) with 
U{xi — xi,X2,x^ — X3) which corresponds to a shift within the plane (xijXa), for the upper 
half-plane we obtain: 


RU{xi - xi, 0 +,X 3 - X3) ■ cr 2 (xi, 0 +,X 3 ) 


.J(x2=f)+) 

—RE2{xi - xijO+jXg - 3:3) • ii(a;i,0+,a;3) 


dxidx3 


J(X2=0+) 0x2 

+/3^ / <P{x[-xi,0~^,x^-X3)q^^\xi,0~^,X3) dxidx3 

J{X2=0+) 

’^3 - —xixi, 0 '^,X 3 ) dxidX 3 

J(X2=0+) 0 x 2 

-F/ 3 + / <^(3:3 - 3:1,0+,X3 - X3)j^^^(a;i,0+,X3) dxidx3 = 0 . 

J{X2=0+) 


(43) 

^ 2 : 2 = 0 +) 

An analogous expression is derived for the lower half-plane where 0+ is replaced with 0“. Sub¬ 
tracting this expression from the (|43p . we obtain 


/ 1 l^[[U'!iix'i - Xi,X3 - X3) ■ (cr 2 )(xi,X 3 ) 

4(x2=0) 


-f 

-f 


R{U){xi - xi, 3:3 - X3) ■ [[0-2]](3:1,3:3) - R{E2){x\ - xi,X3 - X3) ■ [[ii]]( 3 :i,X 3 )| dxidx3 


'(X2=0) 


dd>. 


'yt 


dx2 


(xy - 3:1, X3 - X3) [[6»]] (3:1,3:3) -h 


d^. 


'yt 


dx2 


(xi - Xi,X3 - 3:3)(6 I)(xi,X3) 


+ [[^I 3 t]]ix'i - Xi,X3 -X3) (^2) (3:1,3:3) + (^/ 3 t)(x[ - Xi,X3 - X 3 )[[g 2 ]](a;i,a: 3 )} dxi 

+ i 


{X2=0) 


(xi - Xi, 3:3 - X 3 )[[x]]( 3 :i,X 3 ) -F 


(xi - Xi,X3 - 3:3 )(x)(3:i,X3) 


+ P/3c]](a;i - a:i,3:3 - X3) (j2) (xi, 3:3) (^^c)(a;i - a:i,Xg - X3)[[j2]](®i,^3)} dxidx 3 = 0 , 

( 44 ) 

where the notation and indicates the following normalized potentials: 

<^7c = 7 c^, d>f3t = / 3 t^, ( 45 ) 
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The physical quantitities u,cT2,0,q2,x^j2 involved in the can be represented as follows 


/(X1,X3) = /^+^(X1,X3) + /( ^(X1,X3), 


( 46 ) 


where the superscripts and ^ ^ denote functions whose support is restricted to the positive 
and negative semi-axes, respectively: 


= f{xi,X3)H{xi), /( \xi,X3) = f{xi,X3)H{-Xi), 


(47) 


and H denotes the Heaviside function. The reciprocity identity (j44p then becomes 
i?[[t/]] @ ( 0 - 2 )(+) - r{U2) ® [N^“^ 

= - i^[[^7]]®(o■2)(-)-i^(^7)®[[c^2]](-)-i^(^7)®[[c^2]]W+i^(^2)®M^+^ 




7t 


dxo 


® [[01 


(-) _ 




Jt 


dx2 


® ( 0 ) 


(-) _ 


d<P^ 


7t 


dx2 


® [[01 


(+) _ 


d<P. 


-yt 


dx2 


® ( 0 ) 


(+) 


- ® (^ 2 )^-) - ® [fe]](-) - ® ( 92 )^+) - ® [[<72]]^+) 




7 c 


dx2 


® [[x]] 


(-) _ 


d<P 


7 c 


dx2 


® (x) 


(-) _ 




dx2 


-)® [[x]]^+^ - 


d<P 


7 c 


dx2 


® (x)^+^ 


- M ® - i^^c) ® [[J2]](-) - M ® (J2)(+^ - ® [b2]](+t 


(48) 


where the symbol ® denotes t he co nvolution with respect to both xi and X 3 , which is defined 
as follows (jArfken and Weberl . [2003) ^ 


f ®9 = 


+ 00 


f{xi - XI,X3 - X3)g{xi,X3)dxidx3 = 0 , 


(49) 


On the basis of representation ()46p . ( 0 - 2 )^^'^ is the traction along the interface, ahead of the 
crack tip, and [[w]]^“^ is the crack opening (displacement discontinuity across the crack faces). 
[[[/]] and (U) are the syi nmetrical and skew- sym metrical weight functions matrices defined and 
derived in closed form in Piccolroaz et al.l ( 2009l l. whereas the term {^ 2 ) stands for the traction 
along the Xi—axis corresponding to the singular auxiliary displacements U. 

The integral eq uation (1481) is the general i zation to t hermoelastic diffusive media of the Betti 
identity derived in Piccolroaz and MishurisI ( 2013 1 and Morini et al. (2013a), and it relates the 
physical solution u,cr2,0,x to the auxiliary singular solution U,S 2 - Note that (ilHl) is valid 
for the most general case of static interfacial crack problems in presence of thermal and diffu¬ 
sive effects, and includes also the case of imperfect mechanical interface conditions, associated 
with the discontinuity of tractions and displacements at the interface ahead of the crack tip 
(IVellender et al.l . l2013ll . In the sequel, perfect contact conditions at the interface for xi > 0 will 


be assumed for mechanical fields, then: [['u]]l+l = [[cr 2 ]]^'’'^ = 0 , whereas according to dehnitions 
reported in Section 13.11 discontinuity of the thermodiffusive quantities along the whole plane 
X 2 = 0 will be considered (imperfect thermodiffusive interface conditions). Moreover, the nota¬ 
tion (|26p is used to indicate the symmetrical and skew-symmetrical mechanical loads applied at 
the crack faces, then ( 0 - 2 )*'“^ = (p) and [[cr 2 ]]^~^ = [[p]], respectively. The integral identity 
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becomes: 


RIU]] ® (cr2)(+) - R{U2) 
d<P^t 


d^-yt 

dx2 

d^yc 

dx2 


® [[(9]] - 


® [[xE - 


dx2 

d<P. 


7 c 


dx2 


g [[u]](-) = -i?[[C/]] @ (p) - R{U) ® [[p]] 
® (^) “ {^l3t) ® [[^lE “ [[^/3t]] ® ( 92 ) 

® ix) - {^Pc) ® [[J 2 E - P/3cE ® {k) , 


( 50 ) 


where the terms (0), (x), ( 52 ) and (^ 2 ) stand respectively for the mean values of the temperature, 
mass concentration, heat flux and particle current on the plane 0:2 = 0, dehned respectively by 
expressions (IHUjl and p3]l . Similarly, [[0]], [[xE, [[psE and [[j 2 E denotes the jumps of these functions 
across the same plane X 2 = 0, and defined by (IHTI) and 

The Betti’s identity (1501) will be used further in the text in order to formulate the illustrated 
crack problem at the interface between dissimilar thermodiffusive materials in terms of singular 
integral equations. Explicit integral identities relating the loading applied at the crack faces, the 
temperature, the mass concentration, the heat and mass fluxes to the resulting crack opening 
and traction ahead of the crack tip are derived for the two-dimensional case. 


4 Integral identities 


In this Section, starting from the general identity (1501) . the problem of a two-dimensional inter¬ 
facial crack in presence of thermodiffusion is reduced to a system of explicit integral equations. 
The solution of these equations provides exact expressions for the crack opening and for the 
tractions ahead of the tip corresponding to an arbitrary loading configuration and to the tem¬ 
perature and concentration profiles of the system, obtained by solving the Poisson’s equations 
©• Since thermal conduction and diffusion are both assumed to be isotropic, only normal 
stresses are induced by temperature changes and mass diffusion (see constitutive relation (fT|)). 
For this reason, in the considered two-dimensional problem antiplane stress and deformations 
are not affected by thermodiffusion, so that only the plane strain case is studied in order to 
estimate the contribution of thermal and diffusive stresses on interface fracture phenomena. 

For plane strain deformations, the Betti identity (|5U|) relating the physical solution u = 
[u\,U 2 Y' and <T 2 = [(T 2 i,cr 22 ]^ with the weight function U and U 2 becomes 


i?[[c/E * (cr2)(+) - r{S2) * = -R[[U]] * (p) - R{U) * M 


d^. 


yt 


dx2 

d<Pyc 

dxo 


* [[^E - 

* [[xE - 


d<P. 


yt 


8x2 

8^. 


8x 


* (0) - (^/3i) * [[^lE - * (q'2) 

^ ^ * (x) - * [b2E - P/3 cE * (j2) 


(51) 


where (p) = [(pi), (P 2 )]'^) [[pE = [[biEi [[PsEl^j and the symbol * denotes the convolution with 
respect to the variable xy. Here and in the sequel of the article we will used the following 2x2 
matrices: 


R = 


■ -1 0 ■ 


■ 1 0 ■ 

II 

0 

1—‘ 

_ 1 


1 - 

T—1 

0 

, I = 

, F = 

0 1 


0 1 

-1 0 


1 0 


(52) 
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For plane elastic bimaterials, two linearly independent weight functions, = [C/^, , U 2 = 

= 1,2 are needed in order t o define a complete bas is of the space of s ingular so¬ 
lutions of the homogeneous problem (see Piccolroaz et ah ( 20091 ) and Morini et al. ( 2013bl l for 
details). The weight functions tensors may be constructed by ordering the components of each 
independent solution in columns: 


U = 


Ul 

u} 


ul 


= 


^2 


uf 


(53) 


In order to express the weight function matrix (I53l) (i^ in terms of elastic Lame potentials 
), it is important to note that in the considered plane strain case the vector potential iF 
possesses only one non-zero component directed along 0 : 3 —axis, and then: W = 1^363 = 1 ^ 63 . 
Two couples of elastic potentials ((^^,lF^) and (^^,<F^) are introduced, and the components of 
U are expressed as follows: 


U = 


+ 


m = 




Ut = 


+ 


Ul = 


q^2 


(54) 


dxi dx2 ’ dx2 dxi ’ dxi 8x2 ’ 8x2 8x1' 

Let us introduce the Fourier transform of a generic function / with respect to the variable 


xi: 


/ OO -j POO 

/(xi)e*«"idxi, fix,) = = TT / 

-OO J — OO 


(55) 


Applying 


to expression d^T]) and considering representations (f54)l . we get 

[[UfR{a2)+ - {IJ2)^R[[ur = -[[UfRiP) - (U)^ 

8^ ,1 _ _ 

«) - M (i> 

s2lT 



(56) 


where ^ 7 c = [7c^^ ^/3c = [/3c^\whereas 

the superscripts "*■ and “ denote “ -|- ” and “ — ” functions respectively. In order to obtain from 
expression (|56p a system of integral equations relating the crack opening and the traction ahead 
of the tip with applied loading and thermodiffusive quantities, explicit expressions for the elastic 
potentials and are needed. 


4.1 Derivation of elastic potentials 


Applying the Fourier transform to equations (I54p . the following system of non-homogeneous 
ODEs is derived: 


( §' - = Ul 

\ +i^§ = U2, 


(57) 


where ' denotes the total derivative with respect to X 2 - Since both the pairs of transformed 
elastic potentials (^^,<F^) and (^^,lF^) can be obtained by solving the system (1571) considering 
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respectively U = and U = [^^1,^/1]^ non-homogeneous parts, the superscripts ^ 

and ^ has been omitted. 

Th e Fourier transform of the singular displacements Ui and U 2 was derived in 
( 20n9l V Considering the upper half-plane xi > 0, they are given by 


Piccolroaz et al. 




U2{^,X2) 


X2 - 


—I 


A+ -b 2/r+ 
l?l(A+ +/U+ 

sign(^)x2 


^21 + * 






C(A+ + /i+)_ 


_<e(A+ + /i+) 
T '21 — X2 T 


- sign(^)x 2 

A+ -b 2 /i+ 
|^|(A+ + /r^ 


a;, 


22 


a;. 


22 


2n+ ’ 

= -|Ck2 

2/i+ 


(58) 

(59) 


where ^2ii0 ~ ‘^2i(?!0 ) ~ ^2i(?)0''') "^ 22(0 ~ T'or,(£,0 ) = T'oo(£, 0+) are Fourier 


(2013. 


»)■ 


For the lower 


transform of the traction at the interface (see for example iMorini et al. 
half-plane, replacing —1^| with |^|, //■*■ with and A"*" with A“, similar expressions are found. 
The singular displacements (15811 and (IM|) are substituted into (l57|) . and the solution of the 


linear system yields the following expressions for the transformed elastic potentials, valid for the 
upper-half plane (see Appendix B for details): 


m,X2) 

^{i,X2) 


^ + 2/i+) 


1 

2^ 


X2^J-'^ 

le|(A++ //+)_ 


A, 


'22 


+ H+) 


■U 21 + i 


sign(0 g;2(A+ -b 2/r+) ' 


e-ISp2 

~2^'' 




(60) 


Similarly to the case of the displacements, replacing —1^| with |^|, with and A^ with A“ 
in the (|60p the elastic potentials for the lower half-plane are obtained. 

It is possible now to derive the jump and the average of the Fourier transforms of the 
normalized elastic potentials (I45p across the plane containing the crack. The traces of the 
expressions (I60p on the plane 0:2 = 0 containing the crack are given by 


^(^,0+) 


r 1 1 

0 

1- 

r ^ 

>1 

1 _ 

[ 4^2J 

. ^22 _ 


, !^(C,0+) 


0 , 


i sign(g) 



1 - 

1 


. ^22 _ 


(61) 


^(^,0-) 



a:. 

a:! 


21 

22 j 




0, 


^sign(g) 

4^V 



■^21 


. -^22 _ 


(62) 


so that, taking into account the two linearly independent weight functions defined by means of 
equations (|53p . we derive the following matrix form for the traces of and 


■ <^ne,o+) ■ 


^2(^^0+) 



S. 


1- 


■ ^H^,o-) ■ 


_ ^2(^,0-) . 



A 7 . 


1- 1 


21 ^22 

^2— 

^21 ^22 


yl- V-.1 

^21 ^22 

y2— y-i2— 

^21 ^22 


nl- 



= (i;2) V(0, 


= (a;2) V(0- 


(63) 


(64) 
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Using equations (1631) and (16411 together with the definition (I45p the jump and the average of the 
Fourier transform of the normalized potentials become: 

[[^f,t]] = - Piv-) = {^2f[[ml (65) 

+ I3rv~) = ( 66 ) 

= (i;2)^(/3+r?+ - p-rj-) = {IJ2nVcl (67) 

= i(i;2)^(/3+77+ + /3-r/-) = {E2f{Vc), (68) 

where: 


M = 


1 


0 

ht 


iVt) = ^ 



0 

he 




(69) 


and the expressions for the bimaterial thermodiffusive parameters ht,it,hc and ic are reported 
in Appendix A. 

In order to derive explicit integral identities by equation (l56l) . also expressions for the jump 
and the average of the Fourier transform of the heat and mass fluxes across of the plane containing 
the crack are required. The derivative respect to X 2 of the transformed elastic potentials in the 
upper half-plane (IGOD is given by 


^'{^,X2) 


■x 2 /i+sign(^) /!+ 

L 

X 2 il^ 

A+ 3 / 1 + 


e-ICh 2 

A+-F/t+ ^(A+-F/i+)_ 

^21 + 

_A+ /!+ 

2 |^|(A+ + /!+)_ 

^22 1 

2 / 1 + 


^\C,X2) 


/ 

■x 2 (A+ -F 2 / 1 +) 

A+ -h 2 / 1 + 

T~ i 

■a:2(A+ -F 2 /i+)sign (0 

A+ -h 3 / 1 + 

1 

A+ -h /!+ 

iei(A+ + /x+)J 

2^21 * 

A+ /!+ 

2 e(A+ + /!+)_ 


replacing —1^| with |^|, with fj.~ and A+ with A“ in these expressions, the derivatives with 
respect to X 2 of the trasformed elastic potentials for the lower half-plane are obtained. The 
traces of (p and iF on the plane X 2 = 0 then become 


^'(^,0+) = 

i 

A+ 3/1+ 

^21 

2(A+ + /i+)C’ 

4/i+(A+ + /i+)|^| J 

- ^22 - 

lF'(e,0+) = 

A+ -h 2/1+ 

i(A+-h3/i+) 

r ^21 

2/i+(A+ + /i+)|^r 

4/1+(A+ -h /i+)C . 

1 ^22 - 

^'(C,o-) = 

i 

A“ -|- 3/1“ 

^21 

2(A“ -|- fJ-~)C 

4/i-(A- -F/i-)|^| _ 

[ ^722 J 

lF'(C,0-) = 

A" -h 2/1" 

i(A -|- 3/1 ) 

^21 

- 2/i“(A“ -1- /!“) ^1’ 

4/i-(A- -h/i-)^ - 

- ^22 - 


and the following matrix representation can be introduced: 


( 70 ) 

( 71 ) 

( 72 ) 

( 73 ) 


■ ^i'(e,o+) ■ 


^2'(^^0+) 



U. 


1 - 

21 

2 - 

21 


j^l- 1 
■^22 
y2- 

^22 J 


i 


2 (A+ +^+)C 
A+ -h 3 / 1 + 


(i; 2 )^C+( 0 , ( 74 ) 


4 /i+(A+ + /i+)|C| J 




























































(75) 




^21 

yl- 

^22 

^21 

y2- 

^22 


i 


2(A +/i 
A“ + 3/i“ 

. 4/x-(A- + ^-)|C| . 




The jump and the average of the Fourier transform of the derivatives of the elastic potentials 
across the plane X 2 can be hnally written as 


dx2 


{S2f{-ltC-ltC) = {^2V[[Ctl 


(76) 


(9^, 


it 


dx 2 / " 2^^'^ WC+liC) = {^ 2 Y (C), 


(77) 



(i;2)^(7+C+-7e-C”) = (i^2)^[[C]], 


(78) 


a# 


dx 


\ = -(i;2)^(7+c+ +7c"C-) = (i^2)^(Cc 


(79) 


where: 


[[CJ] 

[[Cc]] 


1 


imt 

rit sign(^) 


1 


irric 

Uc sign(^) 


(C) 

(Cc) 


1 

1 


ipt 

9tsign(0 J ’ 

iPc 

gcsign(0 J ’ 


(80) 

(81) 


and the bimaterial thermodiffusive parameters mt,nt,pt,qt,mc,nc,Pc and Qc are defined in Ap¬ 
pendix A. 

The obtained expressions for the jump and the average of the Fourier transform of elastic 
potentials and fluxes across of the plane 0:2 = 0 are now used together with the Betti identity 
(j56p for reducing the interface crack problem to a system of explicit singular integral equations. 


4.2 Explicit integral identities 

Substituting expressions (I65|) - (l68]i and (IToP - dT^ in equation (l56|) . we obtain: 

gjfR{&^)+ - = -[[c/fB® - 

- {i;2>’’«,)[[91 - (W)^I<,11{9) - {i:2)^(>7,)B2ll - {Sifm («) 

- (S2>^(C.>1R11 - {'h) ■ (82) 

Multiplying both sides by we get 

{,T 2 >+ - BW = -® - AIIpI - (9,)M - fol{9) - (d,)!®]] - [[d.]] ( 92 ) 

-{SojKll - [[9J<X) - (dojCsl - M ih) , (83) 
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where A an d B are the fo l lowing matrices, identical to those introduced in Piccolroaz and MishurisI 
( 2 OI 3 I I and Morini et al. ( 2013al l for the elastic problem without thermodiffusion: 


A = R-^[[Ur^(UfR, B = R-^[[Ur^(B2fR, 
the vectors (g/-), [[firj], (dt) and [[dt]\, associate to thermal effects, are given: 

{g,)=M{Q, M=M[[CJ], {dt)=M{rj^), [[d*]] = M[[r,J, 


(84) 


(85) 


and hnally the terms (gc), [[fl'cE) i^c) and [[dc]], corresponding to the contribution of the mass 
diffusion, are dehned as follows: 


(5,) = M(0, [[9J=M[[CJ, (d,) = M(r?,), [[d,]] = M[[g,l 

where the matrix M is given by 


M = R-^IU]]-^{B 


( 86 ) 


(87) 


The matrices A, B and M can be easily compu ted using the symmetric and skew-symmetric 
weight functions derived by Piccolroaz et al. ( 2009l l : 

m = - idsign{^)E]{IJ2), {U) = - i7sign{C)E]{B2), (88) 

where a and d/b are the so-called Dundurs parameters, b and 7 are bimaterial constants dehned 
in Appendix A. Using expressions (f 88 ]l together with the dehnitions ([ 8 ^ and (l87|l . we get 


A = 


2(62 - d2) 


[{ab — 'yd)! -|- i{ad — 76 ) sign(^)Tl], 


B = - 


lei 


52 — d 2 


[ 6 / -h idsign(^)Tl], 


M = - 


lei 


62 -d2 


[bR + id sign (^)F]. 


(89) 

(90) 


Substituting (f 6 ^ and (f87)l in equation (f85]) . the explicit expressions for the vectors [[d^]], (dt), ^g^]] 
and (g^) are obtained: 


[[dt]] = - 


ht 


(62-d2)|^| 


id sign(^) 

6 


[M = 


1 


i{dnt — bmt) sign(^) 
bnt — dmt 


b‘^-d? 

Similarly, for [[dc]], {d^}, [[gfcE and (g^) we get 

h. 


, (dt) = - 

, iOt) = 


2 ( 62 -d 2 )|^| 
1 


id sign(^) 

6 


2(62 _ ^ 2 ) 


i{dqt - bqt) sig-n.{i) 
bqt - dpt 


[[dc]] = - 


(62-d2)|^| [ 


id sign(^) 

Id) — 

id sign (,^) 

6 

’ ^ 2 ( 62 -d 2 )|e| 

6 


[[ffcE = 


1 


62 -d 2 


i{dnc — bnic) sign(^) 

'a ) — ^ 

i{dqc - bpc) sign(^) 

buc — drric 

’ 2 ( 62 -d 2 ) 

bqc - dpc 


(91) 

(92) 

(93) 

(94) 
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Since explicit expressions for all the terms of the identity (18311 have been derived, we can 
now apply the inverse Fourier transform to this equation, obtaining two distinct relationships 
corresponding to the two cases xi < 0 and xi > 0: 

iP)+T-\,[Am]+T-\,[{9m]]+ 

+-^.■<0 [M (®)]+j^-'<o[(9c)[[x]]] 

+-^.■<0 [idMh]]] + J--'<0 [KE (h)] = (95) 


[{dt)M2]]] + J^-io [[[dt]] m + -F-io[(9 


-1 




+ -^.-Vo[W(x)] 


a;i>0 


wm] + [Ki {*)] = 


(96) 


Note that the term {cf 2 )~^ in the identity (l83|) cancels from the (IMI) because it is a “ + ” function, 
while {p) cancels from the ()96p because it is a “ — ” function. To proceed further, the following 
inverse Fourier transforms are evaluated by means o f the convolutions theorem and distributions 
theory ( R oosl . 1969; Brvchkov and Prudnikov . 1989h : 


Ml?l/(0] = — * -F ^[sign(0/(0] =- — * fixi), 

TTXl OXl TTXl 


(97) 


r 


-1 




m 


= --signxi * f{xi), T ^ 


Ki) 


1 


=-( 7 e«z + In |xi|) */(xi), (98) 


TT 


where the function / can be both the jump and the average of 0,x,Q2 and j 2 , and 'yeui is the 
Euler’s gamma constant. Using the inverse transforms (I97p - (|98p in equations (j95p and (j96p we 
finally obtain the explicit integral identities for plane cracks problems at the interface between 
dissimilar thermodiffusive elastic materials: 


(p>++e<,r'[[«ii+eU’(«)+r>sr’fei]+( 92 ) 
+eSr’M+eir’(x)+I’Sr’fcii+i-U’ m 


xi < 0, (99) 


(^2)<+> + + elf’M + sTV) + r’Sr’M + I>ir’ (®) 

+el?wi'->+ gt*\x) +(92) = (m 


; _F(M_) —)■ F( R ) and : .F(R_) —)• F(R 4 -) are respectively singu l ar ari d 

compact matrix operators ( Gakhov and Cherski . 19781 : Gohberg and Krein . 1958 : Krein . 1958l i. 
where F{M±) is some functional space of functions dehned on R±. These operators are given in 
the form: 




1 


2(62 - d^) 


[{ba - d-i)I + {da - h-i)ES^% = - ^2 _ ^2 “ dE], (101) 
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( 102 ) 




b{da — 67 ) 


S(c) ^ _ 


62-^2 




2(62 - d2) 

The sing ular operator = V^SV- an d th e compact op e rator = V+SV- are defined in 
details in Piccolroaz and Mishurid ( 2013l l and Morini et al. ( 2013al l by introducing the singular 
operator of the Cauchy type S\ 

‘5/ = * /, (103) 


and the orthogonal projectors V±-. 

V±f = 


TTXl 


f{xi) ± Xl > 0, 
0 otherwise. 


(104) 


The vector operators 


^(s+),(s~) 

5 ^IC 


and Q 


(s+),(s-) 
2 c 5 


associate to the contri¬ 


butions of the temperature and of the concentration on the plane containing the crack X 2 = 0 , 
depend by 5^®“) and 


G 


is-) 

It 


G 


is-) 

Ic 


G 


is + ) 
It 


G 


is + ) 

Ic 


1 

{dqt - 6pt)5(® ) 

nis-) 

2(62 _ ^2) 

bqt - dpt 

! ^2t 

1 

[dq, - bpc)S^^-^ 1 

gis-) 

2(62 - h2) 

bqc - dpc 


1 

r {dqt - 6pt)5(*+) 1 

gis+) 

2(62 _ ^2) 

bqt - dpt 


1 

{dqc - 6pe)5(*+) 1 

gis+) 

2(62 - h2) 

bqc - dpc 

! =^2c 


1 

{dnt — bmt)S^^ 

b‘^-d‘^ 

brit — drrit 

1 

[duc — bmc)S^^~^ 

0 ^ 

to 

1 

to 

buc — dnic 

1 

{dnt — 6mt)5(^+) 

b‘^-d‘^ 

bnt — dmt 

1 

{dvtic — 6nc)5(^+) 

0 ^ 

to 

1 

to 

buc — dnic 


where 5^® )= and = P+5 are singular integral operators. \ ^ 

and X> 2 c 5 related to the heat and mass fluxes at the interface are given by 


(105) 

(106) 

(107) 

(108) 
^Ic 


is-) 

it 

■ ■ 

q^is ) _ 

! ^2t — 

ht 

■ ■ 

It — 

1 

to 

1 

to 

6J(^-) 

2(62 - d?) 

bj(^-'> 

is-) _ 

ic 

■ ■ 

! ^2c — 

he 

■ ■ 

Ic 

to 

1 

to 

6J("-) 

2(62 - d2) 

bj(^-^ _ 

is + ) 

it 

■ d/C("+) ■ 

'r)is+) _ 

! ^2t — 

ht 

■ ■ 

It — 

1 

to 

1 

to 


2(62 - d?) 

65^("+) 

is + ) 

ic 

■ ■ 

'r)is+) _ 

! ^2c — 

he 

■ d/c("+) ■ 

Ic — 

4(62 - d?) 

6j("+) 

2(62 - d2) 

bJ^^+'> _ 


where JC and are the integral operators 

ICf = sign(xi) * f, Jf = --{-/eui + In |xi|) * /, 

TT 


(109) 

( 110 ) 
( 111 ) 
( 112 ) 

( 113 ) 


20 

































































and = P_/C, = V+K, = V^J, = V+J. 

Equations (1991) and (jlOOp form a system of integral equations relating the crack opening 
and the traction ahead of the tip to the mechanical loading applied at the crack faces and to 
the values of temperature, mass concentration, heat and mass fluxes on the plane containing 
the fracture. The solution of (j99p by the inversion of the matrix operator provides crack 
opening that corresponds to arbitrary loading configurations acting on the faces and arbitrary 
profiles of the thermodiffusive quantities. The obtained expressions for can then be used 

in (jlOOp for evaluating the traction ahead of the crack tip. 

Note that, in the case where the Dundurs parameter d vanishes, the operators (I1U1I) - (I102|) 
can be written as 


and then the integral equations (fM]l can be decoupled and reduced to: 

- j'Si'iSgld = (p) + _ 2 e5(.)) iip] + + e(--) 

+T>v (®)+eSr’M+sir’fa)+I’ir'fci+ -pt' (*>. 


m+pV 

xi < 0 . 


(114) 


[[92]] 

(115) 


In the case of homogeneous material, we additionaly have a = 0 for the operators (11141) . qt = 
mt = Qc = tUc = 0 for the (|105p - (ll08l) and ht = he = 0 for the (I109p - (I112I) . and the integral 
equations ()115l) become 


_T(.)£!!!2lT 

b dxi 

_V)2!!!2lT 

b dxi 


= (pi) - 


- §S‘- 


Pc. 


= (re) + ^'*>l|Pll + jW + y{x) - 


Xi < 0 . 


(116) 


Assuming symmetrical mechanical load ([[pi]] = [[P 2 ]] = 0) together wit h neg li gible thermodiffu- 
sive effects, equations (I116p become identical to expressions derived by Bice ( 1968l l. 

Differently from the general case (I99p , where the i nversion of the matrix operator by 
means of the general procedure reported in Veknal ( 1971)1 ') is re quired, for the solution of equat ions 
(jll5p and (jll6l) only the inversion of the operator 5^*^ ('see IPiccolroaz and MishurTi (|2013l ll is 
needed. In order to allow some simple examples of application of the obtained integral identities, 
in the next Section uncoupled integral identities (11151) will be solved for certain illustrative cases, 
and effects of temperature and heat fluxes profile on crack opening and traction ahead of the 
tip will be discussed in details. 


5 Illustrative examples 

In this section, the derived general integral identities are applied in order to study the effects 
of two particular temperature and heat flux profiles on crack opening and traction ahead of the 
tip. The temperature and heat flux functions used in the illustrative examples here proposed are 
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derived in Appendix D by the solution of the Laplace’s equation in a semi-plane considering two 
different sets of boundary conditions. Note that, due to the same form of the associate vector 
integral operators (see expressions (Iin5p - (lin8p and (lin9l) - (|112p L the results here reported for the 
temperature field can be immediately generalized for studying the effects of mass concentration 
and mass flux variations at the interface. 


5.1 Punctually localized heat flux profile at the interface 

The following punctually localized prohles for the avegare and the jump of normal heat flux 
across the plane X 2 = 0, containing both the crack and the interface, are assumed: 


{«> (x-i) = - h ) 


hiVxi)^^l,kt + K) 


xi — ai 

I 

xi + ai 

L 


-(5 


-5 


Xi - 02 

L 

Xi - 02 

I 


(117) 

(118) 


where oi, 02 , L > 0 and oi > 02 . It can be easily verified that the functions (11171) and (I118p 
satisfy the integral balance conditions introduced in Section [21 The average and the jump 
of the temperature at 0:2 = 0 corresponding to the heat flux profiles (I117p and (|118p can be 
evaluated solving the steady state heat conduction equation for both the upper and the lower 
thermodiffusive half-space. This procedure is reported in details in Appendix 1C.11 and the 
resulting average and jump of the temperature are given by 


(9) (xi) = [ln(xi - aif - ln(xi - 02 )^] , [[^IKxi) = 0. 


(119) 


In Fig. [21 the variation of the normalized average temperature ([119p and of the normalized 
average heat flux (111711 is plotted as a function of the spatial coordinate xi/L assuming ai/L = 4 
and 02 /^ = 2. 

For simplicity, the Dundurs parameter d involved in the weight functions (1881) is assumed to 
be zero, and then the integral identities (IMP and (HOOP decouple and equations (IMP take the 
form (ES]). Considering only the contributions due to the fluxes and temperature profiles, the 
uncoupled integral identities (I115P become 


5(s) 

Si^) 


dxi 

(9x1 


xi < 0 

-nt(0) + ^j(-)(g2) + 


( 120 ) 


Remembering the definitions of 5^*“^ and introduced in the previous Section, and applying 

these operators to the heat and temperature prohles (I117P , (I118p and ([119p , the following results 
are obtained: 

5(^-)(6() = |ln [(xi - fli)^] - In [(xi - 02 )^] } = 0, (121) 
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Figure 2: (a): Normalized temperature profile at the interface X 2 = 0 reported for ai/L = 4 and a 2 /L = 2; (b): 
Normalized average heat flux profile at the interfaces X 2 = 0 reported for ai/L = 4 and a 2 lL = 2. 


0, 


xi — ai 

z 


-5 


xi - a2 

z 




( 122 ) 


where Tt = lt{kf + k'j~) + ht{kf — k^), and the explicit derivation of expression (jl2ip is reported 
in Appendix D. 

Substituting the (|121|) and (|122l) into equations (|12Up we get 


— = 0 , 


dxi 


5 W 


a[[u2]]^ ) Os 


Xl <0 


(123) 


dxi 






The si ngul ar integral ope r ator is now inverted following the procedure illustrated in IPiccolroaz and Mishuris 
( 2013l i and iMorini et al. ( 2013al l. This leads to 


dxi 


dxi 




+ arctan 


-^ — arctan 

Xl 



(124) 

(125) 


which after integration gives explicit expressions for the crack opening associate to the heat 
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( 126 ) 


fluxes (I117p and ()118p and to the temperature profiles at the interface (I119p : 


* \xi'j — —f - - 


TT 


nt 


TT 


-(oi - 02) + {y/-aiXi - y/-a2Xi) 


ai 


+{xi — ai) arctan w-(xi — 02 ) arctan w- 


Xl 


a2 


Xl 


(127) 


The tractions ahead of the crack tip are given by the solution of equations (|1UUI) . which in 
this case assume the form: 


0 oxi b 

1 nt , . k 


Xl > 0 


(128) 


(a22)(+) = 




Applying the operators and to the flux and temperature functions (jll7p . (jllSp and 

pi9p . the result is given by 

{In [(xi -aif] - In [(xi - 02 )^]} = -y[iJ(ai-xi)-i 7 (a 2 -xi)], (129) 




{s+) 


Xl — fll 

I 


-5 


Xl - 02 

L 


= [(^1 “ ®i)^] “ [(^1 “ “ 2 )^] } , (130) 


where H{ai — xi) and H{a 2 — xi) are Heaviside step functions (see Appendix D for details 
regarding the derivation of (11291) 1. Substituting expressions (I129p and (|13UI) into equations 
(11281) . the tractions ahead of the crack tip hnally become 


(n2i)^^^(xi) =- - H{a2- Xl)], (131) 

(fj 22 )^+^(xi) = —^^nt - {In [(xi - ai)2] _ In [(xi - 02 )^] } . (132) 

In order to study the influence of the different thermal expansion coefficients of the materials 
on the crack opening and on the tractions ahead of the tip, the upper and lower half-space media 
are assumed to have identical values for the elastic moduli A+ = A“ = A, = fi~ = fj,, and for 
the thermal conductivity coefficients = kt, and dissimilar values for 7 ^”. in 

this case, the bimaterial parameters involved in expressions (11271) and (11321) become: 


mt 


It 


2(A -|- /i) 






2(A + /i) 




24 






































where = 7 ^'*'/ 2 (A + /r). Substituting these explicit constants in equations ()127p and ()132h . 


we obtain: 




(133) 


i + \ 

It 


TT 


-(«! - 02 ) + {yj-aixi - y/-a 2 Xi) 


+(xi — oi) arctan ./—— — (xi — 02 ) arctan .— , (134) 

V xi V 3:1 


=-■ * * 


2b 


1 - ^ ) [i7(ai - xi) - H{a2- xi)], 


(135) 




1 + ^ ) {In [(a:i - ai)^] - In [(xi - 02 )^] } . (136) 

7t 


Note that, as expected, for the case of homogeneous material, corresponding to = 7 ^”, 
(<X 2 i)^^^ vanishes indipendently of the results derived by applying the integral operators to the 
temperature profiles. This means that for the case of a plane crack in an homogeneous elastic 
themodiffusive material, no shear stresses are due to the temperature and heat flux profiles on 
the crack surface. 

The variation of the normalized crack opening (I134p is reported in Fig. [3] as a function 
of the spatial coordinate xi/L assuming a single value of the ratio % I= 0.8 and different 
localization of the Dirac Delta heat flux functions (uni) and (111811 . corresponding to different 
values of the normalized distances ai/L and 02 /^. The values for these distances have been 
chosen such that oi > 02 , as imposed in the definition of the heat flux profiles (|117p and (|118p . 
The crack opening profiles shown in Fig. [3)/(a) have been computed considering the same value 
of oi/L = 0.8 and four different values of 02 /^ = {0.001,0.1,1,2}. It can be observed that, as 
a 2 /L decreases, and then the contribution to the heat flux functions depending by 5((xi — 02 ) /L) 
approaches the crack tip, the crack opening increases. In the author’s opinion, this is due to 
the fact that, for small values of 02 /T, the peak of the heat flux profiles (|1I7I1 and (|1I8I1 and 
the maximum of the temperature (IllOp are localized near to the crack tip. Consequently the 
maximum of the thermal load provided by these distributions of heat flux and temperature 
approaches the crack tip, the crack process is favored and an increase of the crack opening is 
pro duced, in analo gy to what is observed for mechanical loading functions localized near to the 
tip ( Freund . 1998lb 

In Fig. [3]/(6) the crack opening [[rt 2 ]]^~^ is reported for the same value of 02 /^ = 2 and four 
different values of ai/L = {2.1,2.5, 3,4}. Differently from what is detected for the variation of 
02 /.b, it can be noted that the crack opening becomes larger as the normalized distance a\/L 
increases. As we expect, [[ 02 ]]^”^ is almost zero for ai/L = 2 . 1 , since in this case oi ~ 02 , and 
then observing expressions (11171) . (11181) and (|119p it can be easily deduced that {qi) ^ 0, [[gi]] « 0 
and {9) ^ 0. Looking at the structure of prohles (I117p . ()118p and ()119p . it can be easily deduced 
that if the value of 02 /^ is fixed and a\/L increase, the effort of the temperature and heat fluxes 
becomes more relevant. Consequently, the loading action provided by the temperature and the 
flux at the interface increases and gives rise to higher values of the crack opening. 
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The variation of the normalized traction (cr22)*''*'^ is plotted in Fig. S] as a function of the 
spatial coordinate xi/L for the same single value of the ratio = 0.8 and the same sets 

of values of ai/L and 02/T assumed for the crack opening. In order to compute profiles of 
the traction ((T22)^'*~\ the symbolic computation program Mathematica was used. As it can be 
expected, the tractions exhibit a singular behaviour at the crack tip, for xi = 0 , and for xi = ai 
and xi = 02, where the temperature profile (jllOp diverges. 




Figure 3 : (a): Variation of the normalized crack opening [[u2]]*'~^ with the spatial coordinate x-ijL, reported for 
different values of a\IL\ (b): Variation of the normalized crack opening [[w2]]^~^ with the spatial coordinate xijL, 
reported for different values of a^/L. 




Xi jL Xj jL 

Figure 4 : (a): Variation of the normalized traction {022)''^'^ with the spatial coordinate xijL, reported for 
different values of a\IL\ (b): Variation of the normalized traction (cr22)*'^^ with the spatial coordinate x\/L, 
reported for different values of a2/L. 

In Fig. [S] the crack opening [[tt2]]*'“^ and the traction ahead of the tip are reported as 
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functions of xxjL for single values of ai/L = 4 and of 02 /-/! = 2 and several valnes of the ratio 
— {0'8j 1-2,1,0.8}. It can be noted that greater values of "it correspond to a greater 
crack opening. This means that, taking into account the considered heat flnx, the temperature 
profiles and the geometry of the model, the crack process is favored if the thermal expansion 
coefficient of the material occupying the lower half-space is greater with respect to that of the 
medium occupying the npper half-plane. 




Fignre 5; (a): Variation of the normalized crack opening [[^2]]^”^ with the spatial coordinate xtjL, reported for 
different values of the ratio "ft j"it i (b): Variation of the normalized traction (0-22)^^^ with the spatial coordinate 
xijL^ reported for different values of the ratio "it j"it ■ 


The traction expressions ()136p can now be used f or cal cnlating the stress intensity factors. 
Applying the standard definition reported in Freundl ( 1998ll . we obtain 


Ki = lim = 

XI —>-0 

Kjj = lim y/27rxi{a2i)^~^'’(xi) = 0. 

xi —>-0 




1 + ^ 

It 


(137) 

(138) 


In Fig. E)/(a) the variation of the normalized stress intensity factor Kj with the ratio 02 /T is 
reported for ai/L = 4 and several values the ratio "ftHt — 1-2; 1,0.8}. It can be observed 

that Kj = 0 for a2/L = aijL = 4. This is due to the fact that for 02 /T = ai/L both the heat 
flux and temperature profiles vanish (see expressions (jll7p . (I118p and (I119M . and consequently 
no thermal load is present at the interface and no crack opening is generated. Moreover, it can be 
noted that as 02 /^ decreases, Kj increases monotonically. This behavionr is in agreement with 
what is detected analysing the crack opening, and also in this case it can be explained observing 
that for small values of 02 /T the heat flux profiles (11171) and (11181) and the temperature (11191) 
possess a peak localized at the crack tip. As a consequence, the maximnm of the provided 
thermal load approaches the crack tip and the crack process is favored. In Fig. El/(^)) is 
plotted as a function of ai/L assuming 02 /^ = 2 and the same values "ft /it — 1-2,1,0.8}. 

According to what is observed in Fig. EF(a)) Kj = 0 for a2/L = ai/L = 2. The reported 
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curves show that the stress intensity factor becomes larger as ai/L increases. This behaviour is 
due to the characteristics of the heat flux and temperature expressions (I117p . (jllSp and ()119l) : 
looking at the structure of these profiles, it can be easily verified that if the value of a 2 /L 
is fixed andai/L, the effort of temperature and heat fluxes becomes more relevant, then the 
loading action provided by thermal effects increases. Consequently, the crack process is favored 
and the value of Kj becomes larger as does the crack opening (see Fig. [3)/(6)). Observing both 
Fig. EF(a) and Fig. EF(^), it can be noted that the stress intensity factor increases as the ratio 
'ItHt- This conhrms the behaviour detected observing the crack opening, and it means that 
considering the heat flux and temperature profiles given by expressions (I117p . (|118p and (|119p . 
the crack process is favored if the thermal expansion coefficient of the material occupying the 
lower half-plane is greater with respect to that of the medium occupying the upper half-plane. 

Considering the case of a homogeneous material, where we have 7 ^ = 7 ^ = 7 t, the crack 
opening (I134p assume the form 


- « 2 ) + W-aiXi - \/-a2Xi) 

TT[X + fl) 


ai 


+{xi — ai) arctan , -(xi — 02 ) arctan , - 


Xl 


0.2 


Xl 


(139) 


(140) 


and then the tractions ahead of the tip (I136p become 


(o-2l)^'^^(®l) = 0, 

(<. 22 )<+>(x.) = [(I, - „,)=] - In [(X. - 027 } 


where 


l g[[n2]](-) 

b dxi 


mtOs 
7r(A -I- 2fi) 



- - + arctan 

Xl 


-^ — arctan 

Xl 



(141) 

(142) 


(143) 


Expressions (I14UD . (|142l) and (I143D can be easily obtained deriving the Green’s function cor¬ 
responding to a thermoelastic semi-plane with the punctually localized fluxes ()117p and ()118h 
applied on the boundary. 


5.2 Symmetrically distributed temperature profile at the interface 

The following distribution profile for the avegare temperature across the plane X 2 = 0, containing 
both the crack and the interface, is considered: 

(0)(xi) = ^e-^. (144) 

The temperature is assumed to be continous at the interface, then [[0]](a:i) = 0 for —00 < xi < 
-|-oo. The average and the jump of the normal heat flux q 2 at the interface corresponding to the 
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Figure 6: (a): Variation of the normalized stress intensity factor Kj with the ratio 02/!/, reported for ai/L = 2 
and assuming different values of !It i (t>): Variation of the normalized stress intensity factor Kj with the ratio 
Oi/L, reported for 02/!/ = 2 and assuming different values of the ratio Ht ■ 


average temperature profile (114411 and to the assumption [[0]](xi) = 0 has been evaluated solving 
the steady state heat conduction equation for both the upper and the lower thermodiffusive 
half-space. This procedure is reported details in Appendix 1C.21 and the resulting average and 
jump of the of the noramal heat flux are given by 


(92) (xi) 


-{kt - kt ) 


L 



T^Erfi 




(145) 




—2(/c+ + ) 


L 



■I^Erh 



2xi 

Ly/n 


(146) 


It can be easily verified that the normal heat flux profiles (I145p and (I146D satisfy the integral 
balance conditions introduced in Section [3l The variation of the normalized average temperature 
(jl44p and of the normalized average heat flux (|145p is plotted as a function of the spatial 
coordinate xi/L in Eig. [71 

Also for this example, the Dundurs parameter d involved in the weight functions ([8811 is 
assumed to be zero, and then considering the contributions due to the fluxes and temperature 
profiles, the uncoupled integral identities have the form (I12UI) . Applying the integral operator 
to the average and jump of the normal heat flux (|145p and (I146p . the following result is 
obtained 

^ W = (147) 
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Figure 7; (a): Normalized temperature profile at the interface X2 = 0 ; (b): Normalized average heat flux profile 
at the interfaces X2 = 0 . 


Substituting expression (jl47jl into equations (jl2np . we get 


5<..M = 


5 W 


dxi 

(9x1 




Xl ^1 

■e~^ 


= -9s nt - 


V L 


Tt\ Xl _d 


Xl < 0 , 


(148) 


L 


-e 


w here Tt is the same bimat e rial p aram eter (defined in Sectio n l5.ll Using the procedure reported 


in 


Piccolroaz and MishurisI ( 201,81 ) and Morini et al. ( 2ni,8al ). the singular integral operator 


can be inverted. This leads to 

dxi 

dxi 


I5W-1 


= -9s nt - 


Tt 


Xl 

—e ^ 

1j 


(149) 

(150) 


The results provided by the application of the operators 5*-^^ ^ and ^ to the function 
_ 

{xi/L) ■ e IT are calculated by means of the symbolic computation program Mathematica. The 
crack opening [[ui]]^~^ and [[^ 2 ]]^“^ are derived by integrating respectively expressions (114911 and 
([I50l). 

Applying the operators to the flux functions (|145|) and (I146D . the following result is 

derived 


2 

1 

TI 


2L 


(151) 
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Substituting expression (jl5ip into equations (jl28p . the tractions ahead of the crack tip become 


( 0 - 21 ) 

{(y22)^^\xi) 


mtOs 


Si^+)( —e-^\ —e~^\ 


V L 




V L 


r. 


-rpT'i - TT 


—e“i? — ( —e~i^ ) 


xi _ 
L 


(152) 

(153) 


Also for the example described in this Section, the upper and lower half-space media are 
assumed to have identical values for the elastic moduli A"*" = A“ = A, = /r, and for the 

thermal conductivity coefficients kf = = kt, and dissimilar values for 74 : 7 ^ ^ li~■ In this 
case, the equations (I149p and (I15UI) become 


= e.s^(i - 2 l ') 5 (.)-. 5 (.-)(^-# 
dxi V 7, / Vi 




9 

ys-^t 


dxi \ 7 ^' 

and the expressions (I152p and (I153p assume the form 


i+A)5'*>-‘(ye-a). 


( 0 - 2 l)^’^^(vCl) = 


I . It 


{cr22)^~^\xi) = - 


it 


5("+) 


L 


L 


(154) 

(155) 


(156) 

(157) 


where the quantity is the same defined in Section 15.1 



Figure 8 : (a): Variation of the normalized crack opening [[rti]]^~^ with the spatial coordinate xijL, reported 
for different values of the ratio 'yll'y^', (b): Variation of the normalized crack opening [[112]]^”^ with the spatial 
coordinate xi/L, reported for different values of the ratio !It ■ 
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Figure 9: (a): Variation of the normalized traction (ct 2 i)*'^^ with the spatial coordinate xijL, reported for 
different values of the ratio (b): Variation of the normalized crack traction {(722)^'^'' with the spatial 

coordinate x\IL, reported for different values of the ratio 


The normalized profiles of the crack opening components and [[^ 2 ]]^”^ obtained 

by integrating expressions (11541) and (I155p . are reported in Fig. [ 8 ] as function of x\jL. Four 
different values of the ratio /^i = { 0 . 8 , 1 . 2 , 1 , 0 . 8 } were considered for the computations. 
For — 1) corresponding to the case of an homogeneous material, we have = 0. 

This behaviour can be directly deduced observing expression (11541) . Observing Fig. [ 8 )/(a), it 
can be noted that the crack opening component changes sign at a determinate distance 

x\lL from the crack tip. As it is shown in the figure, this distance is independent of the ratio 
It ht^ it then varies with the value of a/L and h/L. For > 1, [[^* 1 ]]^”^ is positive near 

to the crack tip, and then becomes negative for xi/L < x\/L. This means that, if 7 ^” > 7 t*~) in 
the xi—direction, associated with the propagation Mode II, the crack opens for xxjL > x{/L 
and closes for xi/L < x\lL. The opposite behaviour is shown for 7 f”/ 7 }*’ < 1, corresponding 
to a crack opening for xi/L < x\/L and to a crack closing for xi/L > x}/L, near to the tip. 
Conversely, the component [[^ 2 ]]^”^ illustrated in Fig. [ 8 ]/( 6 ) is positive everywhere for any value 
of ■ It follows that in the X 2 —direction, corresponding to the propagation Mode I we have 

crack opening for any value of the thermal expansion coefficients of the materials. 

The variation of the normalized tractions (cr 2 i)*''*'^ and (cr 22 )^'*'^ is plotted in Fig. [9]as a func¬ 
tion of the spatial coordinate xijL for the same values of Mt assumed for the crack opening. 
In order to compute the integrals involved in equations (I156p and (I157p the symbolic computa¬ 
tion program Mathematica was used. It can be observed that both (< 721 )^'’'^ and {( 722 )^'^'^ exhibit 
a singular behaviour at the crack tip, for xi = 0. This behaviour is associated with the leading 
term of the asymptotic expansion of expressions (jl56p and (|157) for x^/L —»• 0, correspondi ng 
respectively to the thermal stress intensity factors and ( Liu and Kardomateasl . 2005l i. 

Also in this case, exact expressions for the stress intensity factors can be eval uated by mean s 
of the tractions (11561) and (jl57l) . Applying the standard definition reported in Freund ( 1998IL 
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we obtain 


Kj = lim ^/2^TXl{a22)^'^Hxl) 

xi^O 


Kjj = lim ^/2^:xl{a2l)^'^\xl) 

ii —>-0 


STT^/^VIesS^ f 7 ^\ f r(3/4) \ 

46 V vr(l/4)r(7/4)y'’ 

7r'y^VT0sS+ / _ 7 ^\ 

6r(i/4) V 7yy‘ 


(158) 

(159) 


6 Conclusions 

The problem of a quasi-static semi-infinite interfacial crack between two dissimilar thermodiffu- 
sive elastic materials has been formulated in terms of boundary integral equations by means of 
Betti’s reciprocity identity and weight functions. For the case of a plane strain crack, explicit 
integral identities relating the applied mechanical loading, the temperature, mass density, heat 
and mass flux distributions at the interface and the resulting crack opening have been obtained. 
The solution of this system of singular integral equations provides explicit expressions for the 
crack opening and the stress fields associate to any arbitrary harmonic temperature and mass 
density profile and self-balanced heat or mass flux distribution at the interface. 

Illustrative examples of the application of the integral identities to crack problems character¬ 
ized by punctually localized and distributed heat flux profiles on the interface were performed. 
The crack opening and tractions ahead of the tip corresponding to the assumed temperature and 
heat flux distributions were derived by analytical inversion of the singular integral operators. In 
addition, exact expressions for the stress intensity factors were obtained. 

The proposed original method provides a general boundary integral formulation for interfacial 
crack problems in thermodiffusive bimaterials avoiding the use of the Green’s function and the 
related challenging numerical calculations. The integral identities can be applied for studying 
the effects of any harmonic temperature and mass concentration profile as well as of any self- 
balanced heat and mass flux distributions on crack opening and traction fields ahead of the tip. 
This approach can have various relevant applications, expecially in the modelling of fracture 
processes induced by thermal and diffusive stresses at the interface between different components 
of mutlti-layered electrochemical devices such as solid oxide fuel cells and lithium ions batteries. 
Moreover, the derived integral identities have their own value from the mathematical point of 
view, because, to the authors best knowledge, a similar explicit formulation in terms of singular 
integral equations seems to be unknown in the literature. 
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A Bimaterial parameters 


In this Appendix we report explicit expressions for the bimaterial parameters defined in Sections 
2 and 3. The para meters b, d, 7 and a co ntained in the weight functions matrices (I 88 p are the 
same as defined in Piccolroaz et ah ( 2009l i. here expressed in function of \,[i instead of 


b = 


A+ + 2 /r+ 


+ 


A + 2 /i 


d = 


2/r“*'(A~*" + 2 /i (A + /r ) 2(A“*“ + 2(A + /U ) 

//-(A- + /r")(A+ + 2fj,+) - Ai+(A+ + A+)(A- + 2 //") 


a = 


7 = 


// (A + /r )(A“*' + 2 /i+) + /i+(A“*" + A“*")(A + 2 /i ) 
/i+/r“(A“ + fj.~) + /r+/i“(A+ + /r+) 


(160) 

(161) 

(162) 


// (A“*" + 2/i+)(A +/i ) + /i^(A + 2 /i )(A^ + /i^) 

The bimaterial parameters associated with the effects of the temperature and heat diffusion 
are given by: 


hf = 


It 


It 


rrit = 


Pt = 


It 


2(A+ + /i+) 

it 


4A:+/i+ /i- 

_ 7 r 
2(A~ + /i~) 

+ 


h = 


It 


+ 


+ 


It 


4/ct+/i+ 4fcj /i 


^ (A+ + 3/i+) 7 t (A + 3/i ) 

* 4/i+(A++/i+) 4/i-(A-+/i-) ’ 

74 + (A+ + 3/i+) 74 " (A~ + 3/i") 

— 


(163) 

(164) 

(165) 


2(A+ + /i+) ' 2(A-+/i-)’ 4/i+(A+ + /i+) 4/i-(A-+/i-)’ 

and the constant corresponding to the concentration and the mass diffusion presents the following 
form: 

it Ic „ it , Ic 


rrin = 


Pc = 


ADtp+ ADc fl¬ 
it _ Ic 

2(A~ + fi~) 

Ic 

2(A^ + A^+) 2(A“ + /i~) 


4 = 


+ 


ADtfi+ ADcfi 


2(A+ + /i+) 

it 


^ _ 7c+(A+ + 3/i+) , 7c"(A"+3/i-) 

^r. — . , / I I N “r 


Qc = 


4/i+(A+ + /i+) 
7 +(A+ + 3/i+) 


4/i (A + /i ) ’ 
7 -(A- + 3/i-) 


4/i+(A+ + /i+) 4/i-(A-+/i-)' 


(166) 

(167) 

(168) 
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B Elastic potentials in weight functions space 


In this Appendix, the derivation of the explicit expressions for the Fourier transform of the 
elastic potentials (|6njl is explained in details. These expressions are given by the solution of the 
ODEs system (1571) . which is shown here 


r if' - = ill 

\ 3' + = U 2 


(169) 


where considering the upper half-plane xi > 0, the non-homgeneous terms in the system 
are given by the weight functions 


Ui{i,X2) = 


X2 


A+ 2/r+ 


|e|(A+ + F+) 


E 21 + i 




i{\+ + ^l+) 


U2{C,X2) = 


sign(^)x2 + 




e(A+ + F+). 


X2 + 


- sign(^)x2 

A+ + 2/i+ 
|C|(A+ + /x+) 


A7, 


22 


A, 


22 


2 f+ ’ 

3 id ^2 


2 f+ 


The ODEs system (I169p can be rewritten in the matrix form: 

y +Qy = c, 

where y = [3,3]'^, c = [Di, 1 / 2 ]^ and 


Q = 


0 -ii 
0 


The eigenvectors and the eigenvalues of the matrix ()172p are now defined: 

{Q — \I)w = 0. 

The eigeinvalues are easily determined by solving the characteristic equation: 

det(Q - AJ) = A^ - = 0, ^ Ai = |C|,A 2 = -|C|, 

then the eigenvector matrix W = {wi,W 2 )'^ and its inverse are given by 

1 


W = 


zsign(^) -isign(^) 
1 1 




—zsign(^) 1 
fsign(,^) 1 


(170) 

(171) 

(172) 

(173) 

(174) 

(175) 


Applying the standard eigenvalues method and using the matrices (I175p the system (I17ip is 
written in the following equivalent form 


z + W-^QWz = d, 


where z = W ^y, d = W and 


W-^QW = 


1^1 0 
0 -lei 


(176) 


(177) 
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The ()176p consists of two decoupled equations: 


f 4 - l?l^i = <^1 

14 + l?l^2 = d2 


where the non-homogeneous terms d = W are given by 

di = -^[isign(^)f7i - f72], di = ^[isign{()Ui + Li' 2 ]. 
The general solutions of equations (I178p assume the form 


Zl{C,X 2 ) 

Z2{i,X2) 


C7iel«l^2 _ 


ix 2 T’ 2 i + ( sign(^)x 2 + \ 


H. 


22 


g l^l ^2 

1 ^’ 




(178) 


(179) 


(180) 


Assuming that the displacements and the stresses vanish for xi —>■ + 00 , the following condition 
must be imposed: Ci = 0. The expressions for the transformed elastic potentials y = Wz = 
[ 1 ^, in the upper half-plane can then be obtained by means of the following relations: 

!^ = isign(^)(zi - 2 : 2 ), ^ = zi + Z 2 . (181) 

Substituting expressions (I180p in equations pi81l) . we get: 

^{i,X 2 ) = 

^{^,X 2 ) = 


C 2 e 


-I?h 2 


+ -i 


X 2 lJ.^ 


^(A+ + ^+)' 


^21 + 


1 


X 2 y-^ 


ICI(A+ + 


A. 


isign(^)C'2e 151^2 _^J- 


X 2 (A+ -F 2 ;u +) 
|C|(A+ + /r 


+1 


A. 


21 


sign(0 

2^2 


+ 


e-l?F 2 
2 ^+ ’ 
X 2 (A+ -F 2 //+) 
^(X+ + H+) 


U. 


22 


p ICh 2 
2fi+ 


assuming C 2 = 0, the expressions (I60p are finally derived. Similarly to the case of the displace¬ 
ments, replacing —1^| with |^|, with y,~ and A+ with A“ in the (I182p the elastic potentials 
for the lower half-plane are obtained. 


C Temperature and heat flux profiles 

The temperature and heat flux profiles used in Section [5] are now derived by solving the heat 
transmission problem in a semi-plane for different sets of boundary conditions. Considering 
the upper thermodiffusive half-plane xi > 0 , the temperature profile 6 ^{xi,X 2 ) is given by the 
solution of the two-dimensional Laplace’s equation 

A 6 »+ = 0, 2:2 > 0, (183) 

Two different sets of boundary conditions are considered for the solution of equation (I183|) . 
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C.l Punctually localized heat flux profile at the interface 

The following profile for the temperature on the boundary X 2 = O'*' is assumed 


4 -/ 4 -n Oskt 

qt{xi,X2 = 0 +) = 


Xl — Ol 

z 


-5 


Xl - 02 

z 


( 184 ) 


where ai,a 2 ,L > 0 and oi > 02 . It is easy to verify that the flux function (|184ll satisfies the 
self-balance condition ([36]) (i). Substituting expression (11841) into the heat flux definition (j6j)(i) 
and remembering the condition (|35P (ip the following boundary conditions are defined for the 
the temperature function 9~^ 


9~^{xi,X2 = + 00 ) = 0 , 


89+ 


8 x 2 


X2=0+ 


2L 


Xl — Ol 

z 


-<5 


Xl - 02 

z 


(185) 


Applying the Fourier transform with respect to the variable xi, the equation (I183p becomes 

9+" - ^^9+ = 0, (186) 

where ' denotes the derivative with respect to X 2 - The general solution of equation (11861) takes 
the form 

9+{(,X2) = + C2el«l^2^ (187) 

where the constants Ci and C 2 are determined by boundary conditions (jl85p . They are given 
by 

= C2 = 0, (188) 

consequently, the Fourier transform of the temperature profile for this case becomes 


0 ^{C,X 2 ) = 


A. 

2 |C| 




(189) 


(190) 


and the temperature 9+ is obtained applying the Fourier inversion to expression (11891) 
9 +{xi,X2 ) = -^ {In [(xi - ai)*^ + xj] - In [(xi - 02 )^ + xj]} . 

Assuming the same profile (I184p for the normal heat flux on the boundary X 2 = 0“, an 
expression similar to the (11901) is derived for the temperature 9~ in the lower half-plane X 2 < 0. 
The average and the jump of the temperature across the plane X 2 = 0, which are defined 
respectively by expressions (f3nP (q and (f3TP ny are then given by 


9 


(9) (xi) = [ln(xi - ai)^ - ln(xi - 02 )^] , [[^(((xi) = 0, 


(191) 


and the average and jump of the normal heat flux, defined by expressions (l33|)(i) and (fMP py 
finally become 


{qi) (xi) = -^{kt - h ) 
[tei]](a;i) = + k^) 


Xl — ai 

z 

Xl — ai 

z 


-(j 

-6 


Xl - 02 

z 

Xl - 02 

z 


(192) 

(193) 
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C.2 Symmetrically distributed temperature profile at the interface 

The following distributed profile for the temperature on the boundary X 2 = O"*" is assumed 


0+(xi,X2 = 0+) = 


OsXl 

L 


-e 


( 194 ) 


The boundary conditions for the temperature 9 ~^{xi,X 2 ) are given by expression ()194p together 
with the relation ([35]) (i) here reported 

6 '^{x\,X 2 = +oo) = 0 . (195) 

Also in this case, the Fourier transform of the temperature 9 ^[^,X 2 ) is given by the solution 
of equation ()186p . and then it takes the form (I187p . The constants Ci and 6 * 2 , determined by 
means of boundary conditions (|194l) and (|195l) , become 


i^/^9sL‘^ _i 
=-5-se 


C2=0, 


and then the Fourier transform of the temperature profile is given by 

u/n9s^ 


0^{(,X2) = 


(196) 


(197) 


Applying the inverse Fourier transform to (11971) . the following prohle for the temperature in the 
upper half-plane is derived 


9., 


2 2 


9+{xi,X2) = ye ^ 


f 2X1X2\ 
XI cos ( 1 


—xi Im 


2ix-iX2 

e Erfi 


xi -I- iX2 

z 


— X 2 Re 


-|- X 2 sm 

2ixiX2 

Erfi 


. / 2 X 1 X 2 \ 

'\L^) 


Xi -I- iX2 

z 


(198) 


Assuming the same prohle (11941) for the temperature on the boundary X 2 = 0 , the following 
expression is derived for the temperature 9~ in the lower half-plane X 2 < 0: 


9s -"172 
9 (xi,X 2 ) = —e 


f2xiX2\ . f 2 xiX2\ 

xi cos ( I -I- X 2 sm ' ' 


-|-xi Im 


2ix-^ X2 

e Erh 


xi + 1 x 2 

z 


-|- X 2 Re 


2ix-^X2 

e L'^ Erh 


V ) 

xi -I- iX2 


L 


(199) 


The average and the jump of the temperature across the plane X 2 = 0 containing both the crack 
and the interface are then given by 


{9) (xi) = 


9sXi -fi 
—j—e 72 

1j 


[[ 6 »]](xi) = 0 , 


and the average and the jump of the normal heat hux hnally become 


( 92 ) (xi) = -{kt - )y 

M(a^i) = -2ft+ + fcr)x 


(2xj 


V 


- 1 e"7^Erh 


M_iU-f^Erhf^U 


L2 


L 


2xi 

Ly/^ 

2xi 

Ly/7r_ 


( 200 ) 


( 201 ) 


( 202 ) 


It can be easily verihed that hux prohles (I20ip and (I202D satisfy the integral balance conditions. 
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D 


Evaluation of integrals ^ and for logarithmic func¬ 

tions 


In this Appendix the explicit derivation of results (112111 and (11291) is reported in detail. Assuming 
that the temperature profile at the interface is given by expression (I119p and applying the 
operator the following integral expression is obtained 

{In [(xi - aif] - In [(xi - 02 )^] } 

I J-00 -1 xi -1 

Q 

=-^[Iiixi) - I^ixi)], Xi < 0, 01,02 > 0, ai > 02. (203) 

In order to evaluate , a change of variables is operated 

w = t — ai, dw = dt, xi — t = h — w, h = xi — ai, where oi > 0, xi < 0, h < 0, 
and then this term becomes 

Ii {xi) = Ii {h) = V f h < 0. (204) 

J—00 ^ ^ 

Similarly, for calculating the following change of variables is introduced 

s = t — 02 , ds = dt, xi — t = i — s, i = xi — a 2 , where 02 > 0, xi < 0 ^ £ < 0, 


and this lead to 

r+°o Infs^l 

I-{xi) = I^{t)=V -^ds = 7r2, .^< 0. (205) 

J —00 ^ ^ 

As a consequence, substituting the (|204l) and (j205p into equation (j203p . it can be observed that 
the final result coincides with expression (|12ip introduced in Section PS.II 


S^^-\9) = = 0, yh,i < 0. (206) 

Applying the operator 5*'*'*'^ to the temperature profile (I119p . an integral expression analo¬ 
gous to (12031) is derived 

{In [(xi - oi)^] - In [(xi - 02 )^] } 

= -4| r p M(t-a.)-] ,d 

I 7-00 xi -1 xi - f J 

Q 

=- l:^{xi)\, Xi > 0, 01,02 > 0, oi > 02 . (207) 


39 





















In order to calculate if, a change of variable identical to that operated for evaluating is 
defined 


w = t — ai, dw = dt, xi — t = h — w, h = xi — ai, where ai >0, xi > 0. 


It is important to note that, differently from the term , characterized by h < 0, for the 
evaluation of if both the cases h < 0 and h > 0 must be considered. Consequently the result 
is given by 


lf{xf = lf{h)=V 


r*+00 


ln[r(;^] 
h — w 


dw = 


vr^ h < 0 , 
—TT^ h > 0 , 


xi < ai, 
Xl > Ol, 


(208) 


Similarly, for calculating if the same change of variable used for is introduced 

s = t — a 2 , ds = dt, Xl — t = £ — s, i = xi — a 2 , where 02 >0, xi > 0, 
and the result of the integral is given by 

r-+oo 


/ ■ 

-( 


-ds = 


— s 


vr^ £<0, 

— TT^ £ > 0 , 


Xl < 02 , 

Xl > 02, 


(209) 


The sum of the terms (I208p and (|209p becomes 


0 /i < 0, .£ < 0, ^ 0 < Xl < 02 , 

if (h) — if (£) = ^ 271^ h < 0, £ > 0, ^ 02 < Xl < Ol, 

0 /i > 0, £ > 0, ^ Xl > Ol, 

and then substituting this expression into equation (I207p the following final result is obtained 




0 0 < Xl < 02, 

Os 

— -7^ O2 < Xl < Ol, 
0 Xl > Ol, 


( 210 ) 


The equ ation (I210D can be expres sed by introducing the Heaviside functions H(oi — xi) and 
H(o 2 — Xl) I Arfken and Weber . 200RI 1 . and then it becomes identical to the (I129p used in Section 

IQ 


5(^+1 ( 6 ») = -y [H(oi - Xl) - H(o 2 - Xl)], Ol > 02. 


( 211 ) 
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